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ABSTRACT 


Field  patterns  are  drawn  for  a  plasma  filled  waveguide  of  circular 
cross  section  iimnersed  in  an  ?ixial  magnetic  field.  The  plasma  is  asstimed 
to  be  cold  and  collisionless,  and  linearized  field  equations  are  used. 

Besides  the  exact  analysis  the  MHD  approximation  (low  frequency)  and 
the  quasi-static  approximation  (slow  wave)  are  used.  Comparisons  are 
made  among  field  patterns  determined  by  the  various  methods  of  analysis. 
Among  the  field  patterns  at  different  frequencies  and  wave  numbers 
certain  features  are  pointed  out  and  comparisons  are  made.  These  field 
patterns  include  particle  charge  density  and  current  density,  total 
charge  density  and  current  density,  the  electric  field  and  magnetic 
field. 
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FIELD  PATTERNS  FOR  A  PLASMA  FILLED  WAVEGUIDE 

I.  INTRODUCTION 

The  objective  of  this  thesis  has  been  to  draw  field  patterns 
for  a  specific  plasma  filled  waveguide  in  order  that  field  ijatterns 
of  plasma  filled  waveguides  in  general  may  better  be  visualized.  At 
the  beginning  of  Chapter  I  several  assumptions  are  set  forth.  From 
these  assumptions  the  necessary  equations  for  the  fields  are  derived. 
In  Chapter  II,  the  equations  are  organized  in  a  maimer  used  in  arriv- 
ing at  the  field  patterns  and  then  specific  regions  of  the  (  /3  -^  a)  ) 
diagram  are  discussed.  The  (   3  -u)  )  diagram  is  divided  into  th^ 
magnetohydrodynamic  (MHD)  region,  the  cyclotron  resonance  regiqn, 
the  plasma  resonance  region,  the  region  of  cutoff  fields  (propaga- 
tion constant  equals  zero),  and  a  region  with  the  propagation 
constant  intermediate  between  cutoff  and  resonance. 
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II.  EQUATIONS  GOVERNING  FIELDS 
In  this  chapter  the  equations  of  the  fields  are  derived.  Eight  prelim- 
inary assumptions  are  made  concerning  the  plasma  model  and  the  type;  of 
geometry.  The  first  derivation  is  the  dielectric  tensor  which  comes  from 
applying  conservation  of  momentum  principles  to  the  particles.  Under  the 
heading  of  the  determination  of  the  fields  the  transverse  fields  are  found 
in  terms  of  the  axial  fields  and  the  dispersion  relation  is  derived. 
Under  the  heading  for  the  determinantal  equation  boundary  conditions  are 
matched  for  fields  having  no  variation  in  the  angular  direction.  For  this 
case  the  fields  are  written  out  in  explicit  form.  Under  the  heading  of 
cutoff  fields  a  method  of  solution  is  given  for  finding  the  cutoff  frequen- 
cies. A  discussion  of  the  MHD- regime  is  then  made  and  finally  the  theory 
of  the  modified  quasi-statics  approximation  is  developed. 
2.1  Preliminary  Assumptions 

The  field  plots  of  this  thesis  are  for  a  very  restricted  situation; 
therefore  the  development  of  the  equations  will  be  restricted  from  the 
beginning  and  further  restricted  as  the  development  of  the  equations 
progresses.  A  more  general  development  which  this  thesis  parallels 
closely  is  given  in  (l.)  Waves  in  Anisotropic  Plasmas ,  PP.  (l8-2^), 
CH.  9^  ^^   CH.  10.  The  preliminary  assutnptidns  and  restrictipns  are: 
(l)  A  linearized  macroscopic  description  has  been  chosen  for 
the  analysis.  This  means  that  every  time  dependent  field 
quantity  that  is  determined  is  a  small  perturbation  from 
an  unperturbed  state. 
(2. )  The  system,  waveguide  and  plasma  is  homogeneous  in  time 
and  space.  Fields  of  the  form     c^  ^   will  be  assumed 
and  the  plasma  characteristics  do  no  change  with  position. 


6. 


(3.)  The  geometry  will  be  restricted  to  cylindrical  coordinates 
with  the  radius  of  the  waveguide  a=  .05ni. 

(4.)  The  conducting  surfaces  of  the  waveguide  are  assumed  to  be 
perfect  electric  conductors,  i.e.,  n    X  E,=Oatr"a. 

(5.)  The  plasma  description  is  further  simplified  by  assuming  that 

no  collisions  occur  and  that  the  plasma  is  temperate.  A 

temperate  plasma  is  one  which  satisfies  the  inequalities  that 

the  particle  induced  velocity  be  much  less  than  the  random 

thermal  velocity*  which  must  be  much  less  than  the  phase 

velocity  of  the  wave.  The  lowest  phase  velocity  which  is 

considered  is  in  the  magnetohydrodynamic  region  where  the 

phase  velocity  is  as  low  as  I.75  X  10"2  times  the  velocity 

of  light.  For  the  electrons  to  have  a  thermal  velocity 

less  than  one-tenth  the  phase  velocity  of  the  wave  in  the 

MHD  region  the  temperature  of  the  plasma  must  be  less  than 

10^  °K  or  in  terms  of  electron  volts  the  plasma  must  be 

less  than  a  1  ev.  plasma.  All  the  rest  of  the  patterns 

have  a  phase  velocity  greater  than  one-tenth  the  phase 

velocity  of  light. 

*For  a  Maxwell  velocity  distribution  the  most  probable 

velocity  of  a  particle  is  given  by  u,j  =  j  2Kr'  where  k 

'  m 
is  Boltzman's  constant,  m  is  mass,  and  T  is  temperature. 

It  is  by  definition  the  thermal  velocity. 

(6.)  The  axial  unperturbed  magnetic  field  has  a  strength  of  .1023 

webers  /  m^,  * 
(7»)  The  plasma  is  made  up  of  electrons  and  ions,  i.e.,  protons 

which  have  a  mass  of  I.836  X  10^  times  the  mass  of  an  electron. 


7. 


(8.)  The  unperturbed  density  of  each  type  of  particles  is 
1.808  X  10^'^  particles/m-^«  * 

*rhe  values  of  the  quantities  in  assxaDptions(6. )  and  (8.)  were 
chosen  as  such  for  two  reasons.  The  first  is  that  a 
normalization  of  quantities  which  is  used  in  the  derivation 
of  equations  gives  convenient  numbers  to  work  with. 
Secondly,  the  order  of  magnitude  was  chosen  because  they 
are  thought  to  represent  typical  values  for  a  plasma  filled 
waveguide . 
2.2  Determination  of  Dielectric  Tensor 

The  equations  governing  the  fields  are  Maxwell's  equations  and 
a  conservation  of  momentum  equation  for  each  type  of  particles. 

\7  X  E^  =  -J  ^jUo   H.  (1.) 

V  X  H  =  J  uJCo   E,  /i^  (2.) 

JuJfflk  XJc  =  %(l^/vkXBo)     (3.) 
where  the  k  subscript  stands  for  the  type  of  particle. 

Since  there  is  imiformity  in  the  z  direction,  the  fields  will  be 
broken  into  their  transverse  (T)  and  axial  (z)  components  with^  =  Bq  iz 
Solving  equation  (3.)  for  jir^  gives 

^kT  ' 


JUJ  -  ekBn  iz     ^ 
-w^  /  /ek  Be*  *== 


m^ 
mk     _ 

""^z    =  fk      Ez  (5.) 

Jwnk 

In  order  to  have  the  units  stand  out  some  normalization  constants 

will  be  introduced  and  sane  quantities  will  be  defined. 

"^     '    J  (6.) 

%k=  li^  (7.) 
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V 


u 


u 


=  u' 


+  u'^ 


(8.) 


(9.) 


(10.) 


P     "  I»   ■  ^  pi 

where  c  is  the  velocity  of  light  in  free  space  and  a  is  the  radius  of 

the  waveguide,  u^.^^  is  the  normalized  cyclotron  frequency,  Upj^  is  the 
normalized  plasma  frequency,  Up  is  the  effective  normalized  plasma 
frequency,  and  u  is  the  normalized  frequency.  Assumptions  (6.),  (?.), 
and  (8.)  gives  xx^^   =  —3  and  Up  =  i^-. 

Substituting  the  normalization  constants  into  equations  (U.)  and 
(5.)  and  then  writing  the  velocity  in  matrix  form  gives 

—  Ju  -u 


Vk  =_!!L 


°^k  UJn 


ck 


u2  -u^ 


ck 


-J 


2     2 
u   -  u^ 


0 


ck 


u, 


ck 


u2  -u2 


ck 


^J  u 

"^ 2" 

u  -  u 


ck 


-J 


E  (11.) 


In  order  to  determine  J  in  terms  of  E  some  definitions  are 


introduced. 

A 

= 

No     ej,   ^ 

(12.) 

A 

S 

ie     +   xii 

(13.) 

I 

= 

0"  •     E 

(li^.) 

K 

s 

cr 
I    +        ^ 

(15.) 

'^              i  iJJ  r 

a-  is  determined  from  equations  (ll. ),  (l2.),  (13. )  and  (l^J-.). 
From  equations  (12,),  (13.),  and  (lU. ) 

^°  E,   =  N^  (ee^ie  ^     ^±  ;U)     d^- ) 
Substituting  equation  (ll.)  into  equation  (16. )  determines  Of.  The 

OJ   thus  foxind  substituted  into  equation  (15« )  determines  K  the  dielectric 
tensor.  Carrying  out  the  above  operations  give 
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K 


where 


K^_ 

1 

-.  j 

0 

Kx 

0 

1 

"  r 
1 

0 

0 

(17.) 


Kj_=l  — 


f'^   = 


K 


u^p  (u^   -♦-  ^ce  ^ci) 

-J  ^  ^%   C^ce  *  ^ci) 
(u'-u2ee)ru2-u2^i) 


(X8.) 

(19.) 
(20.) 


J!       Waves  in  Anisotropic  Plasmas  has  the  equations  of  this  section  on 
pages  l8  -  23.  Equation  (11.)  is  equation  (2.8)  on  page  I9,  equation 
(17.)  is  equation  (2.11)  on  page  20,  and  equations  (18.  -  20.)  are 
eqixations  (2.21)  on  page  23. 
2.3  Detenmination  of  the  Fields 

Substitution  of  equations  (l^.)  and  (l5.)  into  eqiiation  (2.) 
leaves  two  vector  equations  with  two  vector  unknowns.  Equations  (l.) 
and  (2.)  are  rewritten  along  with  two  auxiliary  equations  found  by  taking 
the  divergence  of  the  two  curl  equations.  These  auxiliary  equations  are 
not  independent  of  the  first  two,  but  nevertheless,  are  useful  relations. 


V  X   E^   =   -jUijUo  E 

V  X       K       =   JCU^o  K  . 
V  •  K  0  E  =  0 

V  •  H  =  0 


£ 


(21.) 
(22.) 
(23.) 
(2if,) 


Again  the  procedure  is  to  split  the  fields  into  transverse  and 
axial  components 0  Since  uniformity  exists  in  the  z-direction,  the  axial 
dependence  of  the  fields  can  be  taken  as    e    •  The  operators  become 
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V    =     ^^   -    y  r^  (25.) 

^   •   ^   =  K^^  +K^  i;XET  +  Kj|  E2  r^       (26.) 

After  splitting  the  fields  into  their  components,  the  trans- 
versely directed  part  of  equations  (21.)  and  (22.)  can  be  solved 
for  the  transverse  fields  in  terms  of  the  axial  fields.  This  is 
matrix  (27.  )•  Using  this  matrix  and  the  remaining  parts  of  eqxxations 
(18  -  2k)   two  coupled  equations  between  E^  and  H^  can  be  foxind  in 
a  straight  forward  manner.  These  are  equations  (36.)  and  (37' )•  The 
above  procedure  is  carried  out  in  more  detail  in  Waves  in  Anisotropic  Plasmas , 
pages  13*+  -  138»  Matrix  (27.)  is  the  same  as  matrix  (9.21)  on  page  135, 
and  equations  (36.)  and  (37o)  are  the  same  as  equations  (9«30)  and 
(9.31)  repsectively  on  page  I36. 


Jt 


where 


p  1 

1 

R 

«      1    s 

-  —    1- 

T    1 

P 

U        1      Q 

-Q    1 

-S 

P     ]     R 

-U    1 



T     T    p 

i     _ 

^t\ 

^T^z 


n  X   V 


T  \ 


r   X  Vt  Hj 


(27.) 


D=   (  r^  +  xx\^f  ^   (u^Kx)' 
a 

r  =   a  ^ 

p=   -r(r2  +  u2Kj.  ) 


D 


R 


Ju3 


Sc 


Q.   =  Tu^Kx 


S=  ju  (P^-^u^K^) 
T  =  jP^  u  Kx 


> 


€0 


(28.) 
(29.) 

(30.) 

(31.) 

(32.) 
(33.) 

(35.) 
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The  two  coupled  equations  are 


V^  H^  +  f  H^  =  g  E^ 


where 


d  =  j  u  r  Kx 


2 


a2K, 


Mo 


s    ^o 


2k  2 


f=   (r^  -t-  u^  K,  +  u^  K/  \ 


K, 


=    -J  u  r  K^  K, 


a^  K, 


■?o^ 


(36.) 
(37.) 

(38.) 

(39.) 

(^0.) 
(^1.) 


>u. 


Equations  (36.)  and  (37.)  can  be  combined  to  give  two 
uncoupled  fourth  order  equations.  The  equation  for  Eg;  is 

2 


[v^  +  (c  +  f)  Vt  +  (cf  -  dg)] 


Ez  =  0 


(^2.) 


Plane  wave  solutions  of  the  f onn  exp  (-Jp^rr)  satisfy  equa- 
tions (36.),  (37.)  and  (42.)  where  p  is  a  complex  function  of 
frequency  only.  The  application  of  the  identity  (43.)  to  equation 
(U2.)  gives  the  dispersion  relation  for  the  medium  equation  {hk,). 


Vrj,  e 


-J  ^.r^  --J^ee"'^  '^  '^r 


(43.) 


p  -  (c  +  f)  p2  +  (cf  -  dg) 


0 


(44.) 


In  order  to  satisfy  the  boundary  conditions  at  the  waveguide  wall 
two  two  dimensional  Helmholtz  equations  will  be  solved  using 
equation  (42.)  to  give  the  relation  between  the  two  P  valves. 
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(Vt^  +  Pi^j  E^i  =  0  (45.) 

Ez  =  ^zt     +  Ez2  (46.) 

Hz  =  hj  Ez|  +  h2  E22  (U7.) 


where  h^  is  determined  from  equation    (36.) 

\  s    c  -  p| 


(U8.) 


d 
The  above  equations  can  be  found  on  pages  ( 15U . )  and 

(155»)  in  Waves  in  Anistropic  Plasmas . 


2,k   DetermimitsiL  Equation  for  Cylindrical  Coordinates 

For  cylindrical  coordinates  equation  (45.)  is 
^2E3i  +  1   ^  E^i   +  1  ^2  Ezi  +  Pi^  Ezi  ^  Q 


Applying  a  separation  of  variables,  eqxiation  (50. )>  "to 
equation  (49.)  gives  equations  (51.)  and  (52.).  mj^  is  the 
separation  constant  and  is  not  to  be  confused  with  the  m 
used  for  mass  earlier.  The  m  for  mass  will  not  be  used  in  any 
of  the  following  discussions, 

Ezi   =  Ai   R(r)  $((15)  ^^^^^ 

d^I 


(Pir)2  d^R(r) 


d4)"'    +    m,$     =0  (51.) 


In  order  to  have  the  boxmdary  conditions  satisfied  for  all 
^      m^  =  m,  and  for  the  fields  to  be  periodic  with  0    m  must 
be  an  integer.  Finally,  R(r)  -  Jin^^r)  since  there  is  no  need 
for  singular  solutions  in  the  waveguide.  Therefore,  the  solution 

13. 


i 


sought  is 


£ 


zi 


=  A,  J^  (P^r)  e^^^ 


(53.) 


For  the  rest  of  this  section  it  is  assumed  that  m  =  0. 
For  m  :^  0  see  page  226. in  Waves  in  Anisotropic  Plasmas .  The 
boundary  condition   E^  =  0  at  r  =  a  gives 

where 

Xi  =  Pi  a  (55.) 

The  boimdary  condition  that  E^=  0  at  r  =  a  gives  the 

determinantal  equation.  E  ^^  is  determined  from  matrix  (27.)  and 

equations  (^6.),  (^7»)  and  (53«)'  The  determinantal  equation  is 

Xi   Jl(Xi)    _   Xg    Ji(X2J 


where 


bi   Jo(Xi)       bg    Jo^Xg) 
bi  =  1  -  /r^  +  u2  K_L 


(56.) 


(57.) 


a  D 
Rewriting  equation  (^^. )  in  its  explicit  form  gives 


A      1 
4  '  ' 


1  (c  +  f)  -   (c  .  f)V4dg 


(58.) 
(59.) 


The  simultaneous  solution  of  equations  (56.),  (58* )  and  (59*)^ 
for  varying  values  of  u  determines  the  (  fi-u))  diagram.  Fig.  (L.)  is 
the  part  of  the  diagram  that  is  complete  for  m  =  0.  This  diagram 
results  from  the  work  of  Messrs.  P.  Serafim  and  P.  Chorney. 

Using  the  following  four  definitions,  the  complex  functions  for 
the  ^    and  H  fields  for  m  =  0  are  given  by  equations  (6U.)  through 
(69e). 


Ik, 


\-     '\'o(^,    i)    e^"(-n*>     e       '  (^(t]] 


(60.) 


(X3    £)  .  ^"^-^^i  .  -tai)]  (a.) 


^J". 


(620 


'2     = 


^F, 


br 


(63.) 


£ 


E 


0 


i=l 
2 

i=l 


Gi 


'(ruK^)  2     +fp2  +  ^2ji^j2     K         Kjr^+u\)Xi^ 


(6U0 


HKx 


•(p2     +     u2k_^)       +(^^^ 


(65.) 


F^      +     ^2 


(66.) 


H 


\2      /n2      .      2 


Xi 


(67o) 


i=l 


fc    O  T- 

1  o      ^ 


i=l 


2  2         2  o  2 

(p        +  u     kJ        +  (u^  K^) 


(680 


H        =       -^£ 
z  vU  N       F 


i 


i=l 


,2    .    .2 


(P      +uX)  K„        -     K_^Xi' 


ju  PKj^ 


(69.) 


H 


By  inspection  of  equations   (65c)  and  (66o )  give 

^     FT, 


o 

>Wo 


JP . 


(TOO 
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A  more  general  statement  of  equation  (TO*)  <^^^  ^e  found  from 
matrix  (27,,).     It  applies  to  field  patterns  for  which  V^   ^zt 
is  parallel  to  Vt  E^g  »  By  definition 

r^  X    %     =    t  (71„) 

"^T^zi  =  F.  (eJ   i;  (72.) 

Matrix  (27o)  then  gives 
T^   o  Et=  %    .\x    y     ^Q  Vj.  Ei,i  +  S  Vt  Hzij 

^'^  (73.) 

2 
%   ^^     -    %^         ^  (t  V^  E,i  +PVt  H^i)  (^^^) 


A-1 


Therefore       ^ 


>l  ^   (Vi  "^         (75.) 


2o5  Cutoff  Fields  (P  =  0  ) 

By  definition  a  cutoff  wave  is  one  whose  propagation 
constant  equals  zeroo  When  cutoff  occurs  the  coefficients 
P,  Q^  T^  d,  and  g  become  zero  and  equations  (36.)  and  (37») 
become  uncoupled o 

Therefore 

\^       =  (Pe^)    =  ^^  \  ill.) 

X.2  =(p^a)2  .  ,2  (V_^Ki)  (790 


^h  

K 
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E  -  Waves ; 


H,  -0 


"ja 


u 


(80.) 
(81.) 

(82.) 


i^  X  Vj  E^  (83.) 


The  boundary  condition  that  E^  =0  at  r  =  a  requires  X 
to  he  a  root  of  the  Jj^  (Xg)  Bessel  function.  The  cutoff  frequencies 
are  given  by 


CO 

H  -  Waves 


u2««  =  x2_  +  u^ 


\   =  Ab  J^  (Xb  -£-)   eJ-*  eJuC^nt) 


E^  =  0 


Hy=   0 

5c  = 


2 


^o  X'  ^ 


^ 


L^^ 


VtH^  +  i^  XV^H^ 


(8U,) 

(85.) 
(860 

(87.) 
(88.) 


For  the  H^cutoffs  the  bovindary  conditions  become  more 
difficult I  therefore  m  =  0  and  m  =  1  will  be  considered 
separately 0  For  m  =  0  Vm  H  is  in  the  It.  direction.  The 
boundary  condition  E^=Oatr  =  a  requires 


X  V^  H^   =  0 


(890 
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Again  the  bovmdary  condition  requires  X^   to  be  a  root  of  the 

Jn  (X^)  Bessel  functions.  The  cutoff  frequencies  are  given  by 
-•-   h 

equation  (79»)  which  expands  into  equation  (90*  )•  Figxires  (2.)  and 
(3.)  are  plots  of  equation  (79')  and  nine  of  the  solution  points  are 
indicated. 

^c^   +  ^l^co   ■»■  ^2  ^o  +  ^3  =  0  (90.) 

where    ai  =  u^ce  +  ^^c  i  *  ^^%    "*"  ^h  ^^■'■•^ 

ag  =  (u|  -  Uce  Uci)^  +  X^^  (u^/  +  4x  +4)  (92.) 

83=  X^u^je  uci   (u^e  u^i  -  ^%)  (93-) 

2.5J.  ^  -  Cutoff  (m  s  1) 

-^-  (joKr)^y^^!)t,  (94,) 

h  ^   ^ 

Ji  (Xu  il 

■^h  a-        -f 
The  boundary  condition  that  E  ^=  0   at  r  =  a  gives 

Xh  Jq  (x^  )  -  1  =  -J  J^  (95.) 

J  (Xh  )  ^-L 

The  simultaneous  solution  of  equation  (79«)  which  is 

repeated  as  equation  (96.)  with  equation  (95* )  gives  the  cutoff 

frequencies  for  m  «  1 


X^ 

^h 


u2  (k.2  ..    ^J  (96.) 


Figures  (2o)  and  (4o)  for  low  frequencies  and  figures  (3,)  and 

(5«)  for  high  frequencies  demonstrate  one  method  for  solving  for  the 

cutoff  frequencies.  Equation  (95.)  is  plotted  on  figures  (U.)  and 

(5«)<>  The  left  hand  side  of  equation  (95»)  is  plotted  versus  X-^   . 
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The  right  hand  Bide  is  plotted  on  the  same  axis  to  the  same  scale  as  the 
left  hand  side  but  versus  u  .  From  these  figures  a  function  of 
(Xjj   vs.   u)  is  generated  hy  selecting  an  ordinant  on  the  graph  and 
reading  the  corresponding  X^   and  u,  then  repeating  with  a  new  choice 
of  ordinant.  Equation  (96*)  has  been  plotted  on  figures  (2.)  and  (3*)* 
Also  plotted  on  these  two  figures  is  the  (Xh    vs    u)  function  generated 
from  figures  (4.)  and  (5»)o  The  dashed  curves  are  for  imaginary  values  of 
X^.  Therefore^  intersections  of  either  solid  lines  or  dashed  lines 
simultaneous  satisfy  equations  (9^«)  aod  (96.)  and  thus  are  solutions 
from  the  H-cutoffs  from  m  =  1. 

For  finding  field  solutions  in  the  cutoff  region  equations  (83*) 
and  (88, )  are  written  in  a  more  explicit  form.  Equation  (83. )  for 
m  =  0  and  m  =  1  respectively  is 


,(X,^u)   .A,     eJ"^%^'   ,x, 


u 


£0 

or 


(97.) 
(98.) 


Hy  =  j  A,  (Xe,  u)  J^  ^X,  £j  r^ 

H^  .  ,A,(X,,u)  e^^  f,  AlfeJ^  t>-  (j,(x.  ^^]  '  Jj^})%\^^^^ 


'e  a 


^«S 


) 


The  explicit  form  of  equation  (88.)  for  m  «  0  and  m  *  1 
respectively  is 

(100.) 


^      (.,,.)-      ^       ..u(  .„t)   ^^1^. 


^-  Ah   (5^  ,u)   J^(X^  I 


o 


r  \r 


•J  Kx 


^r  -  J  ^. 


^ 


(101.) 
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.jf^ofXhi)    ■       "    ''- 


where  the  characteristics  of  (  — j  Kx  /  K  j_  )  is  given  in  Figures 
(1^.)  and  (5.). 

Waves  ifl  Anisotropic  Plflsmas  considers  the  cutoff  fields  on 
pages  158  -  162  and  page  227. 
2.6  Low  Frequency,  MHD  Regime 

By  definition  the  MHD  regime  is  the  limit  that  occurs  when  the 
frequency  is  much  less  than  both  plasma  frequency  and  the  ion 
cyclotron  frequency.  In  this  limit 
K^  —^   1    .    u/ 


^ce  "ci  (103.) 


K^     ^      -J^  Up^  (^e  +^ci) 


2 2" 

^ce  ^ci 


(10^.) 


K..   ^      /u^  »  2 

- 


-  / ]  (105.) 


(c  -  f)  - 

U   d  g 


->  — 


fl  ( 


IH- 

u^K 

^^ 

Kx 

r 

a  K^ 


(1060 


(107.) 


In  the  limit  of  zero  frequency  (  K  >^  constant),  (K^ ^  O), 

and  (K,  V— cxa  ).  a  comparison  of  equations  (I06.)  and  (107.)  shows 

that      \C^-f)^\    ^^    I  ^     ^&\         >   therefore  from  equations  (58.) 
and  (59°)  P-i  ^^^  c   and  Pp   c::::        f.  Only  a  slight  coupling  exists 
between  equations  (36.)  and  (37.).  The  MHD  approximation  considers  the  two 
equations  to  he  uncoupled^- 

20. 


» 


For  the  H-Wave,  equation  (370  when  (g  =  O),  the  cutoff  frequency 
is  much  greater  than  the  ion  cyclotron  frequency,  hence  it  is  not  in 
the  MHD  region. 


E     -     Wave : 

^'y'       ^z     "^     Pe       ^z     =       ^ 

(108.) 

Xe^     =(M)^     -K*    u2    K^^ 

(109.) 

Equation  (109.),  the  dispersion  relation,  is  plotted  in 
figure  (6.).  In  this  figure,  the  slope  (3n/u  )  approximately 
equals  67,  whereas  a  slope  of  57«1  gives  a  phase  velocity  equal 
to  the  Alven  velocity.  The  definitions  of  phase  velocity  and 
Alvel:!  velocity  respectively  are 


^n 


Va    -     I  "^^^   ^ci^  ^  (111.) 

In  the  limit  of  zero  frequency  Q    and    T,  equations  (32.)  and 

(3^.)*  approach  zero;  and  P    and    U,  equations  (30.)  and  (35.)^  are 
given  by 

P        w  -Ha  (112.) 

r  ^  +  ^"K^ 


U 


-j  a  u  K^ 


(113.) 


giving 


1-^3—    i;     X^tHz  ^^'^-^ 
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Equations  (ll^.)  and  (ll5.)  were  used  In  arriving  at  the  field 
plots  for  the  MHD  regime  jrather  than  using  actual  values  for  the 
coefficients .  The  results  seem  questionable  and  will  be  discussed 
in  the  following  chapter. 
2.7  Modified  Quasi  -  Statics  Approximation 

Near  the  resonances  it  is  possible  to  have  waves  with  arbitrarily 
slow  phase  velocities  in  the  assumed  plasma  model.  By  definition  the 
phase  velocity  of  a  wave  is  its  frequency  divided  by  its  propagation 
constant.  Quasi-statics  is  an  approximation  that  takes  advantage 
of  the  slow  phase  velocities.  In  this  approximation  the  velocity  of 
light  is  assumed  to  be  infinite  as  canpared  to  the  phase  velocity. 
From  this  it  follows  that  the  E^  field  is  approximately  curl  free. 
(This  statement  becomes  more  clear  when  one  assumes  plane  wave 
solutions  of  the  fonn  exp.(  —  J  k,  •  ^)  for  Maxwell's  curl  equation 
for  E. 

-  J  k  X  ^   =   J  "^  *^  H  (116.) 


CI  X   E   =    "^p 
k  c 


-f  „  (117.) 

H 


^  '    H 

For  finite  jH  the  right  hand  side  of  equation  (llT.)  approaches 
zero.)  In  Waves  in  Anisotropic  Plasmas  pages  (162-169)  Professor 
Bers  studies  these  resonances  by  taking  limits  of  the  exact  solution. 
At  plasma  resonance  he  shows  that  the  H  field  approaches  zero. 
For  this  case  quasi-statics  is  a  very  good  approximation.  At  cyclotron 
resonance  he  shows  that  {Ejr  *■  "*  O)  and  Hy  remains  finite.  For  this 
case  the  dispersion  relations  differ  in  the  two  approximations.  The 
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limit  of  the  exact  solution  gives  for  the  normalized  p-number 
2       2 

^"^  k7^  "^  2  ^  '^ii  (118.) 

while  that  found  from  quasi-statics  is 

*^"r^^  (119.) 

Before  proceeding  with  the  derivation  it  should  "be  noted 
that  the  field  plots  are  for  a  modified  quasl-statlcs.  That  is 
the  E  field  as  determined  by  quasi-statics  is  used  to  determine 
an  H  field,  hence  the  name  modified  quasi-statics. 

Waves  in  Anisotropic  Plasmas  has  the  derivation  of  the  quasi- 
static  field  solutions  on  pages  (2U0  -  2^8)  and  for  the  modified 
quasi-statics  on  jjages  (253-25^).  Also  equation  (ll8,)  is  the 
same  as  equation  (9.219)  page  (l64.)  in  the  bopk. 

2.7.1  Derivation  of  Field  Equations 

For  the  quasi-static  approximation  equations  (l20.)  and 
(121.)  are  assumed  to  describe  the  fields. 


V  X   E 


(120.) 


V.  £0  K  .  E   =  0         (121O 


•o 


Since  the  field  is  cure  free  a  potential  can  be  (jLefined 
E    =     —  V(j)  (122.) 
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Splitting  the  field  into  its  canponents  as  previous  proceedures 
followed  give  the  following  equations  from  equations  (l21.)  and 
(122.) 

7t  (J)  +  p2(J)  .  0  (123.) 

P2  =  V^    -^  (12U.) 

j&T  '""^T^^  (125.) 

^z     ^   ^(t>  (126.) 

Figtires  (6.)  and  (7,)  are  a  plot  of  equation  (l2U.)»     The 
"boundary  condition  that     ^z  '     0  at     r  =  a     requires  (X     =     pa) 
to  be  a  root  of  the  Bessel  function. 

The    H-  field  is  determined  fran  the    E  —  field  and  Maxwell's 
equations. 

V.    ^  =      0  (128.) 

For  equation  (l27.)  a  particular  solution  can  be  found  which 
satisfies  the  right  hand  side  and  a  hcanogeneous  solution  can  be 
found  which  will  be  added  to  the  particular  solution  in  order  that 
boundary  conditions  can  be  satisfied. 

Again  the  fields  are  split  into  their  components  and  solutions 
of  the  form  exp.(— J;i^»rjp)  are  assumed.  A  few  algebraic  operations 
on  equations  (127.)  and  (128.)  give 


2k. 


Hpz     =       Y3     E^  (130.) 


where        Y      a    -Ju 


1  '  ':!!  -^       ^K^ 


(131.) 


T^  (132.) 

2  r       W       ^1 


y  _  J    U 


Y3   =     ::lji        LS       k^   k„ 


(133.) 

The  h(anogeneoiis  part  is  foimd  from  equations  (l34.)  and 
(135.) 

V  X  ^    =     0  (13^.) 

^  •  A  ■  °  (.35.) 

Therefore  %      «       —  V  (p  (136.) 

v|  ^^  +  )r^v  =0  (137.) 

and  the  total  H  -  field  is 

H  =H    +    H.  (138.) 

Equation  (138.)  is  complete  when  boundary  conditions  are 

specified.  Since  the  boundary's  are  perfect  electric  conductors, 

the  arbitrary  constant  in  the  homogeneous  part  is  found  by 

lettering  Hy  =  0  at  r  =  a. 

After  defining  some  symbols,  the  E  and  H  fields  are 

written  in  explicit  form. 
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A 


(139.) 


a„=   -iT  =   X 


-ZkT 


fioy 

6. 


m  ■   0: 


E  =   A(u.6,,t) 


K=  A(u,/S^,t) 


-^^    ^l(^  |)    t,    +    J,(X  I)     ij 


(140.) 


U41.) 


(142.) 


(143 . ) 


.-f     ..(XX,     ,J 


+  \h  J.(*  ;>  - 


fin      W      '^'^-M     ' 


m  =   1: 


1  =  A(u,0.(}^t) 


^. 


a 


X     -'i(^I) 


r^  ->- J 


(IM.) 


H=  -i-  A(a,a^.«.t)   ^h,  r^    +    h,T<^    + 


^^z] 


where     h^  :   -i 


P 


h     r 


Y,  j,(x  r)  ^  (Y.+jy^)  iiEll  +  cjh^^ 


xr 


(145.) 
(I4fi.) 


^oiM) 


jYg  JUT)-   (Yi-vjYg) 


Ji(x£j 


a 


+  «h 


Il(finC) 


r 


'^^  a 


•^n  3 


^J,(xi)  + 


S^i^'^n-I)] 


(147.) 


(148.) 
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III.  FIELD  PLOTS 

The  section  describing  the  conventions  has  the  equations  which 
describe  the  i)article  motions  listed  and  a  discussion  of  the  two 
conventions  used  in  generating  the  field  patterns.  The  two  conventions 
are  (l.)  polarization  ellipses  at  various  locations  to  describe  the 
field  strength  and  a  vector  to  describe  its  direction  and  (2.)  spacing 
of  lines  to  determine  the  strength  of  the  field  and  direction  of  the 
lines  to  describe  the  direction  of  the  field.  All  the  rest  of  the 
sections  consider  the  plots  for  various  regions  of  the  (B-uj)  diagram. 
Ccanparisons  are  made  between  various  plots  and  the  interesting 
characteristics  of  certain  field  patterns  are  pointed  out. 


3.1  Equations  and  Conventions  Used 

The  objective  of  this  thesis  has  been  to  use  those  equations 
which  were  developed  in  Chapter  I  in  order  to  arrive  at  patterns 
of  the  fields  for  selected  points  of  the  (Q-w)   diagram.  Figure  (l.) 
has  indicated  on  it  those  points  for  which  field  patbfems  have  been 
drawn.  In  order  to  get  the  equations  for  the  field  patteriB  E^  and 
H  were  determined  frcan  those  equations  developed  in  chapter  one  in  the 
various  sections.  The  rest  of  the  field  patterns  were  calculated  by- 
using  the  following  equations.  From  equations  (l4.)  and  (l5«) 


u  I  ^o 


(Xit9.) 


Equations  (8.),   (ll.),  and  (l2.)  give 


4 


2     2  -2—2       '     ° 


2 


-Uck 


^^^r. 


0 


-J  u 


I 1  -J 

!       . 0        I  IT 
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(150.) 


£  =  ^e  +  J.i  (151.) 

^  '  ^,  V.^^  (152.) 


P   =  --r-T-    V  •  J  (153.) 

Taking  the  curl  of  equation  (2.)  and  equation  (l53.)  give 

:/>  «  V  .  €0  E  (15U.) 

The  constants  for  the  plasma  are  given  "by 

P«  ^^     1836.1   ^ 


«ce  =  -3   ^ci   =  ■,3.  (156.) 

Ib36.1 


Equations  (l51.)  and  (l53»)  were  used  to  check  the  correctness 
of  the  answers.  Equation  (151.)  was  used  to  check  that  sum  of  the 
current  densities  as  detennined  frcm  equation  (150.)  actually- 
equalled  the  current  density  as  found  from  equation  (1^9. ). 
Equation  (l530  was  a  second  method  for  calculating  the  total 
charge  density  and  was  used  to  check  the  correctness  of  the  answers 
as  found  by  equation  (154).  This  check  is  equivalent  to  equation 
(23.)    V  .  K^ .  JE^  =  0.  Finally,  since  the  computations  were 
performed  by  slide  rule,  no  more  than  three  significant  figures 
can  be  expected. 

No  i)atterns  were  drawn  for  the  separate  particles  except  for 
frequencies  below  (u  -  l).  This  is  because  nearly  all  the  current 
is  from  movement  of  electrons  at  high  frequencies.  Therefore,  the 
pattern  for  the  electrons  is  the  same  as  the  total  current.  The  pattern 
for  the  ions  is  nearly  the  same  as  that  for  the  E  -  field. 
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2         2 
The  justification  of  this  statement  comes  from  assuming  u   »  u  . 

and  looking  at  equation  (l50.)» 

The  most  noticable  case  for  disagreement  between  the  patterns  for 
E  and  J^  would  be  for  E»ji  to  be  linearly  polarized.  This  would 
produce  elliptical  polarization  for^j^m  with  the  ratio  of  major 
to  minor  axis  equal  to  1.6  X  lo'"^.  Therefore,  it  is  concluded 
that  the  term  (  u  ci  ^     X  E)  has  negligible  affect  on 
the  field  pattern  and  the  current  density  pattern  for  the  ions  is 
approximately  the  same  as  that  for  the  E  field. 

For  the  field  patterns  two  sets  of  conventions  were  used  to 
give  two  sets  of  field  patterns.  The  first  convention  has  ellipses 
of  polarization  at  various  locations  to  show  the  relative  strength 
of  the  fields.  In  each  ellipse  a  line  is  drawn  from  the  center  to 
the  outside  which  indicates  the  strength  and  direction  of "the  field 
at  time  t  =  o  when  the  constant  of  amplitude  for  the  fields  is 
taken  as  a  positive  real  number.  As  time  increases  the  arrow 
indicates  the  direction  in  which  the  vector  rotates. 

The  three  components  of  any  field  are  drawn  to  the  same  scale 
unless  it  is  noted  otherwise  by  (X   Mumber)  on  the  pattern.  The 
particle  current  and  the  total  current  are  drawn  to  the  same  scale. 
Piwally,  the  scale  of  any  particular  (  E-  J  ,  or  H  )  was  chosen  to 
give  convenient  sized  ellipses  of  polarization. 
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Only  two  of  these  views  are  necessary  since  the  third  can  be 
drawn  from  the  first  two.  The  (r  -  (p)  view  is  taken  looking  down 
the  waveguide  axis  (  /  z -direction)  at  z  -  0;   the  (  r  -  z)  view 
is  taken  looking  across  the  waveguide  into  the  minus  (p-direction 
at  4^  =  0.  The{^(t>-z)   view  actually  is  generated  from  the  preceeding 
views.  It  corresponds  to  looking  into  the  minus  r-direction  at 
(J)  =  0  and  moving  the  polarization  ellipse  up  or  down  a  given  dis- 
tance r  depending  whether  the  ellipse  is  at  plus  r  or  minus  r 
respectively  when  ^  =  0. 

The  second  convention  indicates  the  strength  of  the  field 
by  the  spacing  of  lines.  This  convention  has  been  applied  to  all 
of  the  transverse  fields  for  m  =  1  and  to  the  radial  fields 
for  m  a  0.  The  (^-directed  component  is  divergence  free  for 
m  =  0  and  its  fields  are  readily  visualized  from  the  patterns 
of  the  first  convention.  The  z -dependence  of  the  patterns  was  not 
accomplished  due  to  lack  of  time  but  can  be  sketched  easily  from  the 
patterns  of  the  first  convention.  For  an  increase  in  z  the  vector 
rotates  in  the  direction  of  the  arrow.  When   ^  -    ^   the 
vector  has  rotated  by  90°  in  space. 

The  second  conventdon  illustrates  two  features.  It  clearly 
illustrates  the  divergence  of  the  fields,  especially  at  cutoff. 
Laying  the  jjattem  for  the  E- field  or  for  the  J- field  rotated 
by  90°  iii  space  over  the  charge  density  contours  shows  that  where 
the  density  is  the  greatest  the  number  of  lines  starting  or  ending 
per  unit  area  is  also  the  greatest.  The  second  feature  is  that  the  total 
transverse  pattern  rotates  in  the  counterclockwise  direction  for  m  =  1. 
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Before  proceeding  with  the  discussions  of  the  individual  regions  on 
the  {B'US)   diagram,  it  should  be  noted  that  the  assumed  model  of  the  plasma 
may  break  down  at  the  boundary.  Inherent  in  the  model  is  the  assumption 
of  specular  reflection.  Specular  reflection  accounts  for  the  fact  that 
particles  cannot  flow  through  the  walls  of  the  waveguide  as  the  field 
patterns  indicate, 

3.2  Low  Frequency,  MHD  Regime 

For  this  region,  four  sets  of  fields  were  calculated,  two  by 
the  HHD  approximation  and  two  by  the  quasi-static  approximation. 
Such  good  agreement  occurred  between  the  two  approximations  that  only 
two  sets  of  patterns  were  drawn,  one  for  m  =  0  and  one  for  m  -  1. 
The  largest  variation  between  the  two  approximations  was  a  factor  of  two 

for  JT^^T'  ^^^P' 

For  this  region  and  exact  analysis  was  not  possible  because  of 

failure  to  find  a  solution  point  on  the  (B-UJ)  diagram.  A  comparison 
of  the  coefficients  which  were  used  with  those  from  an  exact  analysis 
using  the  same  frequency  and  propagation  constant  will  indicate  the 
degree  of  the  error.  For  an  assumed  u  =  17  x  10~5  and/a^  ■  1.2  X 
10"^  one  of  the  p  numbers j equation  (58«)*  is  very  close  to  c,  equation 
(38.),  but  for  the  exact  analysis  it  is  not  necessarily  a  root  of  the 
Bessel  function  since  there  are  two  p  numbers.  Coefficients  which  would 
have  been  used  for  the  approximate  analysis  to  find  the  transverse  fields 
are 

P  =   j  12,3  (158.) 

Q  -       OoO  (159.) 

T  =   0.0  (160.) 

u  =  1.53  (161.) 
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while  the  coefficients  from  an  exact  analysis  are 

(162.) 

(163.) 
(i6i*.) 
(165.) 

From  a  comparison  of  the  coefficients  and  knowing  that  there 
are  two  p-numbers  each  of  which  are  different  from  the  approximations 
raises  an  element  of  doiibt  as  to  the  value  of  the  approximations. 

Assuming  the  approximations  are  a  valid  representation  of  the 
fields,  it  can  then  be  concluded  that  for  both  m  •  0  and  m  =  1 
the  fields  lie  in  the  transverse  plane  except  for  jJ  and  ^,     The 
strong  axial  electron  current  is  a  result  of  (K^-j — *— -  ©o)  as  {UJ—*-  O). 
Another  interesting  feature  of  this  pattern  is  that  the  ions  move  in  the 
transverse  plane  in  such  a  fashion  that  the  divergence  of  the  total 
current  is  negligible  compared  to  the  divergence  of  either  particle  . 
current,  i.e.,  the  particle  charge  density  is  approximately  3  X  lo3  times 
greater  than  the  total  charge  density. 

3.3  Electron  Cyclotron  Resonance 

Three  sets  of  field  patterns,  two  for  m  s  1  and  one  for 
m  s  0,  were  drawn  for  this  region  frcsn  the  quasi-statics  approxima- 
tion. Here  larger  p-n\mbers  than  those  used  for  other  regions  were 
assiamed  for  two  reasons.  The  first  one  was  to  develop  some  patterns 
with  more  variations  in  the  radial  direction.  The  second  and  most 
important  reason,  is  that  the  dispersion  relation  as  determined  by  quasi- 
statics  approaches  the  true  dispersion  relation  only  for  large  p  numbers. 
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This  can  be  seen  by  a  comparison  of  eqiiations  (ll8.)  and  (ll9.)»  The 
model  for  the  quasi-statics  approximation  assumed  the  velocity  of  light 
to  be  infinite  compared  to  the  phase  velocity  of  the  wave,  in  this  case 
the  phase  velocity  in  the  transverse  direction,  hence  large  p  numbers. 

Two  points  seem  to  be  indicated  from  these  patterns.  The  first  is 
that  some  of  the  lines  for  ^   and  Ht  enclose  upon  themselves.  The 
^  for  other  regions  is  completely  encircling,  cutoff  and  MHD,  or  no 
encircling  lines  as  plasma  resonance.  The  second  point  is  that  a 
ccmparison  of  the  fields  at  electron  cyclotron  resonance  with  those 
at  plasma  resonance  shows  that  the  H-field  for  electron  cyclotron 
resonance  is  about  an  order  of  magnitude  greater  than  the  H  -  field 
at  plasma  resonance  when  the  two  fields  have  similar  propagation 
constants  and  similar  values  for  E^ . 
3.U  Plasma  Resonance 

This  is  an  interesting  region  because  several  comparisons  can 
be  made.  As  a  first  choice  the  exact  analysis  and  the  quasi- 
static  analysis  at  u  s  4.o8  will  be  made.  Here  it  can  be  seen 
that  there  is  very  good  agreement  between  the  two  patterns  and  field 
solutions.  The  difference  which  occurs  is  that  the  E  -  field  becomes 
slightly  elliptically  polarized  in  the  transverse  plane  for  the  exact 
analysis  and  the  strength  of  the  z -directed  field  as  found  by  the 
exact  method  is  about  twice  that  from  the  quasi-static  method. 

The  second  interesting  comparison  is  that  at  similar  propagation 
constants  the  plasma  resonance  cones  closer  to  being  an^-wave  than 
cyclotron  resonance  as  was  discussed  in  the  preceeding  section.  Also 
Eg   for  plasma  resonance  is  about  twice  E^  at  cyclotron  resonance  when, 
the  two  have  similar  magnitudes  for  ^ 
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The  third  interesting  canparison  is  to  compare  the  patterns  of  those 
fields  which  have  solution  points  that  lie  on  the  curve  on  the  {&-\Jj) 
diagram  that  originates  at  ^  Z   3.9U.  See  Fig\are  (l.)«  As  the  propagation 
constant  increases  the  z  component  of  the  E  and  J  fields  grow  while  the 
H^  -  field  decreases.  The  transverse  H  -  field  grows  from  nothing 
at  cutoff  to  strong  radial  then  to  almost  circular  while  the  E  -  field 
starts  from  almost  circular  and  decreases  to  almost  radial.  The^  field 
is  very  similar  to  the  Bt^   field  except  that  it  has  a  strong  radial  field 
at  cutoff. 

The  fourth  interesting  point  is  that  J  and  H  have  a  polarization 
in  the  transverse  plane  which  is  opposite  to  that  of  the  E  -  field. 
The  E  vector  rotates  in  the  counterclockwise  direction  while  J  and  H 
rotate  in  the  clockwise  direction. 
3.5  Cutoff  Fields  (  P  s  0  ) 

.  3.5»1  E  -  Waves 

All  that  seems  apparent  for  these  cutoffs  is  that  the  H  field  is 
the  only  one  with  transverse  components  and  the  charge  density  is  zero. 
They  are  similar  to  the  cutoffs  for  E  -  waves  in  an  empty  waveguide. 
They  differ  in  that  the  cutoff  .frequency  is  higher  than  that  for  the 
empty  waveguide  and  that  the  empty  waveguide  has  no  particles,  hence 
no  current  density.  For  the  empty  waveguide  (^;o-  Xe  )  of  equation 
(8.0)  and  is  a  root  of  the  Bessel  function. 
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3. 5*2   H  -  Waves 

m  =  0: 
Three  field  patterns  were  drawn  for  m  =  0.  For  the  cutoff 
at  u  =   .0U30  the  ion  and  electron  charge  density  nearly  cancel 
each  other.  From  the  patterns  of  ^e^  ^^  iti-T  ^^*  particle  paths 
can  be  drawn  for  low  excitation,  i.e.,  localized  motion  of  particles, 
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(a.)  The  small  ellipse  is  the  path  of  the  ion. 
(b.)  The  large  ellipse  is  the  path  of  the  electron. 

For  u  =  3*9^  "the  interesting  feature  is  that  a  nearly 
circularly  polarized  Et  -  field  produces  a  strong  radially  polarized 
J  -  field.  The  reason  for  this  is  that  Uj.gE^  approximately  equals  u  E^^ 
in  equation  (15O. )  hence  (j^  ^*^0). 

The  cutoff  frequency  u  =  6.3^+  has  nothing  that  particularly 
stands  out.  A  comparison  of  the  three  patterns  shows  that  as  frequency 
is  increased,  the  E-field  goes  from  nearly  radial,  to  almost  circ\ilar 
with  the  stronger  field  being  r-directed  and  then  to  almost  circular 
with  the  strongest  field  being  0  directed.  For  m  =  0  the  H  wave 
at  cutoff  in  an  empty  waveguide  has  the  E  field  entirely  in  the 
0  direction;  therefore  none  of  the  patterns  are  like  those  of  the 
empty  waveguide.  The  higher  the  frequency  the  more  nearly  the  modes 
look  like  empty  waveguide  modes. 
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m  =  1: 

For  m  «  1  there  are  four  sets  of  field  patterns.  A  comparison 
of  the  patterns  for  U  =   ,0U326  and  m  =  1  with   U  »   .0^30  and  m  =  0 
shows  that  for  both  sets  the  charge  density  of  the  electrons  tends  to  cancel 
that  of  the  ions,  that  is,  p^  "  ""  3»75/>  and  p    ■  h.f^p.       The 
particle  paths  for  m  =  1  look  like 


^r 


r  _3       ClIjiLS^^ 


a"°       \      ^      /  a~T  -^  \ 


0 


(a.)  The  smaller  ellipse  is  the  i)ath  of  the  ion, 
(h.)  The  large  ellipse  is  the  path  of  the  electron. 

In  order  to  visualize  the  way  the  electrons  tend  to  cancel 
the  charge  density  of  the  ions  the  second  convention  is  the  better. 
A  comparison  of  u  r  3,37  for  m  »  1  with  u  s  3.9I+ 
for  m  =  0  shows  that  both  have  nearly  circular  polarization 
ellipses  across  the  cross  section  of  the  waveguide,  and  both  have 
the  same  direction  of  polarization. 

The  last  two  sets  of  patterns  are  for  imaginary   p  numbers. 
The  fields  for  U  =  5*07  show  how  large  imaginary   p    nxombers  produce 
very  strong  fields  near  r  =  a.  This  is  because  a  modified  Bessel 
function  grows  rapidly  as  its  argument  grows.  The  fields  for 
u   :  5*77  are  for  p  numbers  that  are  practically  equal  to  zero. 
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3.6  Patterns  Between  Cutoff  and  Resonance 

Two  sets  of  patterns  have  been  categorized  under  this  heading, 
one  for  u  =  U.773  and  one  for  u  »  1.00.  For  the  patterns  at 

u  Z   U,773  one  of  the  normalized  p- numbers  is  8.65  and  the  pther 
is  2.U0.  The  large  p  -number  causes  three  variations  in  the  field 
pattern  as  r  goes  from  zero  to  a.  For  this  pattern  it  is  neces- 
sary to  use  the  patterns  from  the  second  con^vention  because  the  large 
number  of  variations  renders  the  first  set  of  patterns  practically- 
useless. 

The  special  feature  the  patterns  at  u  Z  1  have  is  that  this 
is  the  only  set  which  was  drawn  for  complex  p— numbers.  For  it 
a  new  complex  constant  of  amplitude  was  chosen  in  order  that  the 
fields  would  be  in  either  the  (x"  or  (|>  )-directions  at  time 
t  =  0.  Like  the  cutoff  fields  with  imaginary   p  -  numbers  the 
fields  are  strongest. near  the  outside  of  the  waveguide. 
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MHD: 

E-Wave    (1. )  pp.    169-172 

ra  =  0 

u  =  17  X  10"^ 

S-^  =  l.lij.1  X  10"*^ 

p   =  2,if05/a 

A(z,t)    =  A^eJ^^^xlO-^uint-l.liflxlO-^  |) 

Conplex  Fields : 

E  =  A(2,t)    [-.J795  Ji(pr)t^  +  J^ipr)t^] 

H  =  A(z,t)    [.J103.8  J^(pr)t^] 

J  =  A(2,t)    [23.7  Ji(pr)tp   +  J2.38   Ji(pr)t^ 

Jq  =  A(2,t)    [+1.273x10'"^   J^(pr)t^-j2.25xl0^   Jl(pr)t^ 

-ji4..99xlO"*'3   J^(pr)t^] 
1^  =  A(z,t)   [23.7   Ji(pr}t^   +  j2. 27x10^   Jl(P^)^<j> 

-j 2.  72x10"^  J^{pr)t^] 
p  =  A(z,t)    [-J3. 38x10"'^  J^Cpr)] 
Pq   =  A(z,t)    [-jl.  118x10"^   JqCpx*)] 
Pi   =  A(z,t)    [+jloll8xlO"^   J^Cpr)] 
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^uasl-Staticj   (1.)  pp.  Zl^l-Zk^,    253-251^- 

m  =  0 

u  =  17x10"^ 

/3^  =  5.87x10"^,   /3=  .117[|, 

p  =  2.i4.05/a 

A(z,t)  =  j.^HH^lO'^'^nt-S.QYxlO'^z/a) 

Complex  Fields? 

1  =  A(x,t)    [.jl|10   J^ipr)t^   +  Jo(pr)t^] 

H  =  A(z,t)   [(-6.MpclO"^  J^{pr).3. 35x10"^  |)t^ 

-J103.9   J-L(pr)t^   +j(2.6ipclO"^  J^ipr) 

-1.  ll^xlO"^)t  ] 
z 

J  =  A(z,t)    [12.20  J^(pr)t^   +  jl.230  J^(pr)ti 

-ji^. 99x10^   J^(pr)t^] 
J^  =  A(z,t)    (6.56x10"^   J^(pr)t^  -   j  1.159x10^   J^(pr)t. 

-j4. 99x10*^   J^ipi')t^] 
1^  =  A(z,t)    [12.20  J^(pr)tp   +  jl. 171x10^   J^(pr)t. 

-j2. 272x10 "^  Jjpr)t^] 
P  =  A(z,t)  [-jl.746xlO"7  j^(pr)] 
p^  =  A(z,t)  [-j5.  75x10"^  J^(pr)] 
p^  =  A(z,t)   [+j5.  75x10"^  J^Cpr)] 

Note.'   Field  Patterns  Drawn  For  MHD  Only 
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MHD: 

E-Wave    (1. )  pp.    169-172 
m  =  1 

u   =  17x10"^ 
/3^  =  1.353x10"^ 


P   =  3.832/a 

A(Kz,t)   =  A^eJ(^7xlO-^u;^t.l.353xlO-Va)^j^ 

Complex  Fields: 

E  =  A(<t),z,t)   [j5.93xlO''(JQ(pr)   -       p^      )t^ 

.5.93xl02-l_-t^  Wl(pr)t^] 

J,  (pr)  Jn  (pr)    . 

H  =  A(<>,z,t)   [65.1  -^=p7-t^   +  j65.1(J^(pr)   .-J~^t^] 

Jn  (pr)  ^ 
J  =  A(<|>,z,t)   [-17.63   J^ipr)   +  15.85     \^      )t^ 

J.Cpp)  ^ 
-J  (1.779   J^ipr)  +  15.85  -~F-")^<)> 

-jl^.. 99x10^   J^{pr)t^] 
J^  =  A(<t),z,t)    [-(9.1^9x10-^  J^(pr)   -  1.672x10^  .2-_)t^ 

+jl. 672x10^ (J^(pr)  -     ^p^      )t^ 

-ji^., 99x10^  J^(pr)t^] 

p  Jt  (pr )    . 
J^   =  A((t),z,t)    [(-17.62   J^(pr)   -  1.519x10''  -i_-)t^ 

P  P   J,  (pr ) 

-j(lo697xlO''  JQ(pr)  -  1.519x10^     ^^^      )ij> 

-j2o  72x10"^  J^(pr)t^] 

p   =  A(<t>,z,t)    [«jt.  02x10""^  J3^(pr)] 
Pe   =A((|),z,t)    [-jl.32ipclO"^   J-j^(pr)] 
p^   =  A(<>,z,t)    [+jl. 325x10"^   J-^ipr)] 
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Quasi-Static:      (1.)  pp.2lfl-2l^6,   253-25if 
TO   =  1 

u  =  17x10"^ 

^n  "  9.36x10"^    3=  .1872 
P   =  3.83/a 
A(^>,z,t)   =  A^eJ(17xlO-^u/„t-9.36xlO-3z/a)^j^ 

Conplex  Fields: 

J-,  (pr  )  _.  Jn  (pr ) 

E  =  A(+,z,t)    [ji|.09(JQ(pr)    -       py      )?^  -  k09      \^       t^ 

J-i  (pr )    .  Jt  (pr )    . 

H  =  A(*,z,t)   [65.1  -^p-  t^  +  j65.1(J^(pr)   -  •^=57— )t^ 

+j  1.652x10"^   J,  (pr)t,] 

Ji(pr) 
1  =  A((|),z,t)   [(-12.19   Jo(pr)   +  10.95       py      )t^ 

Jl(pr)  ^ 
-j  (1.228   J^(pr)  +  10.95       ^^      )t^ 

-jl|..99xlo3   J^(pr)t^] 

J«  =  A(^,z,t)   [.(6.5lpclO  -^   Jo(pr)   -  1.156x10'^     \^      )t 

2  Jl^P^>  ^ 

+jl.  156x10^  (J^(pr)  -  -ij— -)t^ 

-ji^.99xlo3   J^(pr)t^] 

Ji  =  A(<j),z,t)   [(-12.19   Jo(pr)  -  1.048x10*^  -J^_.)t^ 

2  .  2  "^i  ^P^ )  ^ 

+^(-1.169x10^^  J^ipp)   +  1.01^8x10''  -i— — )T. 

-2,72x10"^  •^l(P^)^2^ 

p   =  A((|),z,t)    [-j2,77xlO'"7  j^(pr)] 

Pe   =  A(<|),z,t)    [-j9. 16x10"^  Jl(pr)] 

p^  =  A((|>,z,t)   [J9. 16x10"^  Jl(pr)] 

Note:      Field  Patterns  Drawn  For  MHD  Only 
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ELECTRON   CYCLOTRON  RESONANCE: 

Q,uasl -Static  Analysis 

m  =  0 

u  =  2.1+ 

3      =  10.08  ;     ^=  201,6 
n 

P  '■^ 

A(z,t^  =  vJ(2-^^nt-10.08^) 

Coup  lex  Fields: 

i  =  A(z,t)   [-j.547  Ji(pr)t^   +  Jq^P^'^^z^ 

if  =  10"^xA(z,t)   [(1.6i|.l  J^(pr)   +  1.937x10"^!^ (  ,3r))t^ 

-J2.05  Ji(pr)t^  -   j(,898  J^(pr) 

+  1.937x10"^  I^(  6r))t^] 
J  =  A(z,t)   [.3ii4  Ji(pr)?^  -  j.lj-30  Ji(pr)T^ 

-j.35i^  Jo^P^^'^z^ 
p  =  A(z,t)    [-j2. 23x10"*^   J^ipr)] 
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ELECTRON  CYCLOTRON  RESONANCE: 
Quasi-Static  Analysis 
ra  =  1 
u  =  2,k 
Q^   «  7.01  ;  3=  114.0.2 

P  =^ 

A(♦,z^)  =A^eJ(2.WnW.0l|)  ^^ 


Complex  Fields: 

E  =  A(<>,z,t)   [U^kliJjpr)   -  -ip5— )tj,  -  -A?  -^jr-  ^^ 
+  J^(pr)t^] 

H  =  10  ^  A(+,z,t)   [(-2.36  J^ipv)   +  5.31     -^p^ 

-6.ypclO--^(I^(ar)   -     ^^^      )% 

+j(2.95  J^Cpr)  -  5.31  -^ e.iOpclO  ^  -i-j^— )t^ 

-j  (1.292   J^^Cpr)   -  6.ypclO"^  Ii(/3r))t^] 

Ji  (pr )  ^ 
J  =  A(<),z,t)    [(..3iA  Jo(pr)   +  .7iA       pr      )ty 

Ji  (PP)  ^ 
+j  (.14-30  JQ(pr)  -  .7kk       pr      )t^ 

-j.1937  J3^(pr)t2] 

p  =  A(<|),z,t)   [-j  1.611x10"^   Jl(pr)] 
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ELECTRON   CYCLOTRON  RESONANCE; 

Q,ua3i -Static  Analysis 

7.016 


m  =  1  P  =       a 

u  =2,18 

/3n  =  9.93       ^=  198.6 

z 


Conplex  Fields? 

E  =  A((|>,2,t)    [j.707(JQ{pr)   -  -^=57— )t^  -  .707  --?— ^+ 
+  J^(pp)tp] 

H  =  10  ^xA(<|),z,t)   {(-1.1^25  JQ(pr)   +  3.38  -~jr- 
+  1.79ipclO  ^(1^(6 r)   .  -i^_.)t^ 

+j  (1.962  J^(pr)  -  3.38  '—-'  *  l-79ipclO  ^  -4^^— )t^ 
.j(1.003  J^Cpr)  +  1.79ipclO*^  I^{  f^r))t^] 

J  =  A(<|),z,t)   [(-.308  J^Cpr)   +  .732       pp      )t^ 

+j(,i^l^  J^(pr)   -  .732  -^=pF~)^(|> 
-j.390  J^(pr)t^] 
p  =  A((t),z,t)   [-j2.  63x10"*^   Ji(pr)] 
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PLASMA  RESONANCE?      (1.  )pp.    l62-l6if 
Exact  Analysis 
ra  =  0 
u  =  k'Od 

^n  "  ^°-^^ 

P,    =^^ 
*^1         a 

P2  =  i^  =  3^2 


Conplex  Fields: 

E  =  A(2,t)   [-j2.09xlO"^(J^(p^r)  -  e.^BxlO"^!  ^iQ,^r))t^ 
+  3.9ipclO~^(J3^(p3^r)   -  2.4.9x10"^  I^(^r))t^ 
+   (Jq^Pi"*^   "*"  1-^7x10"^  Iq^V^^^z^ 

H  =  A(z,t)  [-2.71xlO"^(J^(p^r)  -  2.i4.7xlO*^  l^{Q^r))t^ 
+J1.722xlO"^(J3^(p^r)  +  3.25x10"^  ^±(^2^^^^^ 
+j6.25xlO"^(J^(p^r)   -  9.89x10"^  ^o^V^^^z^ 

7  =  A(z,t)   [.8.l5xlO'2(J^(p^r)   -  30  36x10""^  I^(Q2r))t^ 
+j5. 17x10 ~^(J^(P3^r)    +  3.9x10"'^  l^{q,^r))t^ 
-j2. 08x10"^ (J^Cp^^r)   +  l.i<27xl0"^  Io^V^^"^z^ 

P  =  A(2,t)   [.jxl0'^(lo9W  Jo^Pl^^   "^  3.76x10"^  Iq((^2^))] 
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plasma  RESONANCE.' 

Quasi-Static  Analysis 
m  =  0 
u  =  i^..08 
a„  =  6,23     ,'      /3=  164.6 

j(1^.08u>^t-8.23  f) 


n 


A(z,t)  =  A^e 

Coup  lex  Fields ; 

E  =  A(z,t)   [-j2. 92x10"^   Ji(pr)"?p  +  J^{pv)t^] 

H  =  A(z,t)    [-5.  69xlO'"^(J^(pr)   -  1.0i|.5xlO"^  \{^r))t^ 

+ji^..20xl0'"^  J^(pr)t. 

+j6.92xlO'^(J   (pr)   -  8.61x10"^  I^(6r))T  ] 
o  o  z 

J  =  A(x,t)   [-1.322x10"^  J^(pr)tp  +  j9. 72x10"^  J^(pr)t+ 

-j.208  J^ipr)t^] 

p  =A(z,t)    [-jloSSxlO""^   J^(pr)] 
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PLASMA  RESONANCE  -  PAST  WAVE: 
Exact  Analysis 
m  =  0 
u  =  i<-.  l^k 
8j^  =  1|..00 
p^  =  3.ll|9/a 
P2  =  jl.87Va  =  j   q^ 

A(z,t)   =  A^eJ  ^         * 

Con^lex  Fields  J 

E  =  A(z,t)   lJ(-.686  J^ipr)   +  .092?  I3^(q2r))t^ 
+   (.393  JiCp^r)  -  .0782  I-^iq2^))t^ 

+  (Jq^Pi^^  ■"  -^^^^  ^0(^2'' ^^^z^ 

H  =  10~^xA(z,t)  [-1.006  J-,^(Pj^r)  +  .1997  I3_(q2r))t^ 
+J  ( .  252  J^  (p-^r )  +  .  0620  I^  ( q2r  )  )t^ 
+J(-795  JqCPi^)  -  .0937  Io(^2^^)^z^ 

J  =  A(2,t)  [-.1717  J-^ip-^r)   +  .OlSkS  l^{q2r))t^ 
+j(.0437  J-^ip^r)   +  .OOij.75  1-^(921- ))t^ 

-j(.2o4  J^Cp^r)  +  .0300  loCqsr))^^.] 

p  =  A(z,t)  [-J10'^{lo090  JqCp^p)  +  ,0728  I^Cqs^))] 
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PLASMA  RESONANCE: 

Quasi-Static  Analysis 
m  =  1 
u  =  ij..26 
Q^  =9.62     ;         8=  192. 1^. 

P   =^ 

^  (1^.26  ujt-9.  62  f) 

Coii5>lex  Fields: 

E  =  A(<>.z,t)   [j.398(J^(pr)   -  -^jr-)^p  -  -398  -^jt"  ^<j, 
+  J^Cpr)!^] 

H  =  10  -^xACKz^t)    [(.i4.97  Jo(pr)   -  .847     \^ 
+  1.129xlO-^[I^(3r)   -  -i_— ])t^ 

-j(.350  J^ipr)   -  .81|.7  -^p7—  -  1.129x10  ^     "^  ^^      )t^ 
+j(.1977  J^Cpr)-  1. 129x10*^  IiCfir))^^] 

J  =  A(<),z,t)   {(.1571  Jo(pr)   -   .268      "^p^      )t^ 

-j(.1108   J^(pi»)   -  .268     "-p^      )t^ 
-J.  1992   J^(pr)t^] 
p   =A(^,z,t)    C-jl. 970x10"^   Jl(pr)] 
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CUTOPP  FIELDS: 
B-Wave 
ra  =  0 
u  =  i^..  67 


-  2A05 

ji^..67iXint 


P  a 


A(t)    =  A^e' 

Complex  Fields: 

E  =  A(t)    [J^ipr)t^] 

H  =  A(t)    [j.1367  J3^(pr)t^3 

J  =  A(t)   [-j.1818  Jq^P^^^z^ 

J^  =  A(t)    [-j.1818  Jo(pr)t^] 

J^  =  A(t)    [-J9. 89x10"^  Jq^P^^^z^ 

p   =0 


E  ^  J 


C-.  (l)=o,z=o  Cr"' 
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CUTOFP  FIELDS: 
E-Wave 
m  =  1 

p  =li82 
*^         a 


A(t,t)   =  A^9^5.5lni*„t  ^J* 

Cod?) lex  Fields: 

1  =  A(<|»,t)   [J,(pr)t J 

i  ^  J,  (pr)  _^  JiCpr) 

5  =  1.833x10-3  A(t,t)    [-  -1^^  t^  -  j(J^(pr)  -  -^57—)^^] 

J  =  A(Kt)   [-J.l53i^  Ji(pr)t2.] 
p  =0 


i 


E  *  J 


t5  _  (1-0,2  0  C'' 
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t^>o,z-o    [^^' 
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CUTOPP  FIELDS; 
H-Wave 
ni  =  0 

p  =i^ 

^  a 

Coup  lex  Fields: 


(1.)     u^^^   =  .01,30        A(t)   =  A^e^0^30vA)nt 

S  =  A(t)   J^ipr)    [2.72x10^  t^  -   jl^..23  t^] 

H  =  A(t)   Jo(pr)t^ 

J  =  A(t)   J^(pr)    [-j.619  t^   +  77.0  t^] 

Jq  =  A(t)   J^(pr)    [+J2.30  t^  -  77.0  t^] 

J^   =  A(t)    J^(pr)    [-J2.92  t^   +  6.56x10"^  ti] 

P  =  A(t)   J^(pr)   [1.83xlO"7] 


(2o)  u^^2  =  3.91,       A(t)   =  A^e-J^-^l^-UJnt 

E  =  A(t)   J3^(pr)    [.1,.98  t^   -   J3.88  t^] 
H  =  A(t)   J^ipr)t^ 

J  =  A(t)    J^^Cpr)    [jlOl^.O  t^  -  l4..1,6  t^] 
p  =  A(t)   J^(pr)    [-3.37x10"''] 


■  •»«a<»«kBaaio*flHioaas9i»a>ca«acS9a«ci«i««aa«B«Jca«i>uiBa«ia«aaa9«4««aims«*i 


(3«)  ^co3   =  ^-3^       A(t)   =  A^e^^-3^^nt 
E  =  A(t)    J^^Cpr)    [310  t^  -   j62i,  t^l 
5  =  A(t)    JoCpr)!^ 

J  =  A(t)   J^^Cpr)   [-jlOlf.6  t^  -  133.2  t^] 
P   =  A(t)   J   (pr)    [2.10xlO"''7] 
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CUTOFF  FIELDS: 
H-Wave 
m  =  1 

u  =  0OI4.326 
A(<>,t)   =  A^eJ-^^33UJnt  3J+ 

Complex  Fields? 

JnCpr)    . 
E  =  A(<>,t)   C-.(272  J^(pr)   -  275       p^      )ip 

+j(if.l7  J^(pr)   -  275-i5p-)t^l 

H  =  A(Kt)    [Ji(pr)^2] 

J-,  (pr )  ^ 
J  =  A(^),t)   [j(.6l5  Jo(pr)   +  76.5       p^.      )ip 

-(76.9   J^Cpr)  -  76.5  -^i=5^r-)^^] 

,    Ji(pr)    . 
J3   =  A(<|),t)    [-j(2.29   Jo(pr)  -  79 . i^-  -=~r- A^ 

.(76.9   Jo(pr)   -  79'k  ~-)t^li 

,    Jn(pr)    . 
Ji   =  A(^>,t)    [j(2.90  J^(pr)   -  2.Qk     \^      )\ 

J,  (pr) 
-(6,52x10-^  J^(pr)    +  2.Qk     \^      )lj>] 
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CUTOPP  FIELDS? 
H-Wave 
ra  =  1 
u  =  3.37 

A{*,t)  =A^eJ3.37u.nt3J* 
Complex  Fields: 


J.  (pr)    . 
E  -  A(*,t)    [5W   J^(pr)   -  23.1  —^^-^^ 

+j(^2i^  J^(pr)   +  23.1  -^=pF->^4^ 

H  =  A((|),t)   J^ipr)t^ 

JtCpi-)^^ 
I  =  A(<>,t)    [-j(98.0  J^Cpr)   -  52.9       pp      )\ 

Hhk'9  Jo^P^^  ""  ^^-^     pF-^^4^ 

p  =  A(<t),t)    [-2.35x10"''^   J-L(pr)] 
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CUTOPP  FIELDS: 
H-Wave 
ra  =  1 
u  =  5.07 
p=j    lk20=j    q 

A(^t)   =  A^e^^-Q7u>nt  ^j^> 

Cornplex  Fields: 

_                                            ^  I,  (qr)  -J  ^ 

E  =  A(<),t)    [(1-985x10^  -i-j: 1.814.8x10^  I^(qr))tp 

-j  (1.985x10-^  -^ 1.366x10^  I^(qr))t^] 

H  =  A(<t),t)   l-^(qr)t^ 

P   =  A((|),t)    [-!{..  57x10"^  Il(qp)] 
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CUTOPP  FIELDS: 
H-Wave 
m  =  1 
u  =  5. 77 
p   =  j.0120ii  .  ^  q 

A(*,t)   =A^eJ^-77UJnt  ^j  + 
Conplex  Fields? 


1  =  51+4  A(<>,t)   [(1   +   (|r)tj.   +  j(l  -   (|)    )t^,] 

H  =  A(<>,t)(|)t^ 

J  =  52.8  A(+,t)   [.j(l  +  3.l6(|)2)t^  +   (1  .  3.l6(|)2)t^] 

p  =  3.85x10"'^  A(<)^t)(|) 
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Discrete  Frequency: 
m  =  0 
u  =  1 

p.  =  IMI  ,f^^ 
c.  a 

,     A(z,t)   =  A^eJ        "  ^ 

CoDplex  Bessel  Functicsti 
JjL(pr)   =  U^(pr)   +  jV^(pr) 

Conplex  Fields: 

E  =  A(z,t)    f  j(-3o55  U^^Cp^r)   +  1.068  Vj^Cp^^r))^^ 

-(.957  U^^Cp^r)   +  1.588  V^Cp^^r))  t^ 

+(l.i^90  U^Cp^^r)   -   .666  V^(p^r))t^] 

if  =  A(2,t)   [(6.22x10""^  ^i(Pi^)    ■*■  10.29x10"^  ^l^^l^^^^r 
-(17.77x10"^  ^l^Pl^^   "^  6.83x10"^  V^(p^r))t. 
=j  (16.27x10'"^  U^(Pir)  +  3.99x10"^  V^(Pi^))^2^ 

J  =  A(2,t)   [(.679  U-^Cp^r)  +  .393  V^(p^r))t^ 

+j  (-1,231  Ui(Pi3?)   +  »172l^  V^(p^r))t^ 
+j  (-1.266  U^(Pir)    +  o566  V^(p^r))t^] 

P   =  A(z,t)    [jxlO""^  (-1.702  U^(p^r)   +  2.26  V^(p^r)] 
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Discrete  Frequency: 
ra  =  0 
u  =  ii-.  773 

5  n  =  1-  ^^3 
p^  =  2.1j.05/a 

Pp  =  8.  65ij-/a 


A(  z,  t)   =  A,e 


j(1^.773uj^t-1.663  f) 


Complex  Fields: 

1  =  A(z,t)   [j(o248   J^Cp^^r)  -  3.15  J^ip^^))"^^ 
+(=.337  Ji(Pir)    +  0643   J^(p2r))t^ 
+(J^(P3^r)   -»•  .666  J^(P2^^y^z^ 


H  =  10"3xA(z,t)    [(.312   J^ 
+j(lo568  JiCp-j^r 

+j(»olj-5l  JqCPiI* 
J  =  A(z,t)    [(0OO96   J-^iPjT 

+3(0 0533   J-^ip^r 
-j{.1779   JqCPiF 


p^r)   -  .B9k  Ji(Pir))t^ 
+  o298   J^(p2r))t^ 

+  3ol9   Jo^P2^^^^z^ 
-  .807  J^(p2r))tj. 

+  .393   J3^(P2r))t^ 

•t-  ,1183   Jo(P2^^)^z^ 


p   =  -^10"°^  A(z,t)    [.1883   Jo(Pir)   +  5.01   Jq(p2I')] 
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